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Abstract 

In  this  paper,  we  address  the  problem  of  the  existence  of  supercon¬ 
vergence  points  of  approximate  solutions,  obtained  from  the  Generalized 
Finite  Element  Method  (GFEM),  of  a  Neumann  elliptic  boundary  value 
problem.  GFEM  is  a  Galerkin  method  that  uses  non-polynomial  shape 
functions,  and  was  developed  in  [4,  5,  24].  In  particular,  we  show  that  the 
superconvergence  points  for  the  gradient  of  the  approximate  are  zeros  of 
certain  systems  of  non-linear  equations  that  do  not  depend  on  the  solution 
of  the  boundary  value  problem.  For  approximate  solutions  with  second 
derivatives,  we  have  also  characterized  the  superconvergence  points  of  the 
second  derivatives  of  the  approximate  solution  as  the  roots  of  certain  sys¬ 
tems  of  non-linear  equations.  We  note  that  it  is  easy  to  construct  smooth 
generalized  finite  element  approximation. 

AMS(MOS)  subject  classifications.  65N30,  65N15,  41A10,  42A10, 
41A30 

Keywords:  Superconvergence,  Generalized  Finite  Element  Method,  Interior 
estimate,  Superapproximation. 


1  Introduction 

The  superconvergence  in  the  finite  element  method  (FEM)  is  a  phenomenon, 
where  the  order  of  convergence  of  the  finite  element  error,  at  certain  special 
points  in  an  element,  is  higher  than  the  order  of  convergence  of  the  maximum  of 
the  finite  element  error  over  that  element.  These  special  points  are  called  natural 
superconvergence  points.  To  the  best  of  our  knowledge,  this  phenomenon  was 

‘Institute  for  Computational  Engineering  and  Sciences,  ACE  6.412,  University  of  Texas 
at  Austin,  Austin,  TX  78712.  This  research  was  partially  supported  by  NSF  Grant  #  DMS- 
0341982  and  ONR  Grant  #  N00014-99- 1-0724. 

1  Department  of  Mathematics,  215  Carnegie,  Syracuse  University,  Syracuse,  NY  13244. 
E-mail  address:  banerjee@syr.edu.  WWW  home  page  URL:  http://bhaskara.syr.edu.  This 
research  was  partially  supported  by  NSF  Grant  #  DMS-0341899. 

^Department  of  Mathematics,  University  of  Maryland,  College  Park,  MD  20742.  E-mail 
address:  jeo@math.umd.edu.  WWW  home  page  URL:  http://www.math. umd.edu/~jeo.  This 
research  was  supported  by  NSF  Grant  #  DMS-0341982. 


1 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

2005 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-2005  to  00-00-2005 


4.  TITLE  AND  SUBTITLE 

Superconvergence  in  the  Generalized  Finite  Element  Method 


6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Office  of  Naval  Research, One  Liberty  Center, 875  North  Randolph  Street 
Suite  1425,  Arlington,  VA, 22203- 1995 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

In  this  paper,  we  address  the  problem  of  the  existence  of  supercon-  vergence  points  of  approximate 
solutions,  obtained  from  the  Generalized  Finite  Element  Method  (GFEM),  of  a  Neumann  elliptic  boundary 
value  problem.  GFEM  is  a  Galerkin  method  that  uses  non-polynomial  shape  functions,  and  was  developed 
in  [4,  5, 24].  In  particular,  we  show  that  the  superconvergence  points  for  the  gradient  of  the  approximate 
are  zeros  of  certain  systems  of  non-linear  equations  that  do  not  depend  on  the  solution  of  the  boundary 
value  problem.  For  approximate  solutions  with  second  derivatives,  we  have  also  characterized  the 
superconvergence  points  of  the  second  derivatives  of  the  approximate  solution  as  the  roots  of  certain  sys¬ 
tems  of  non-linear  equations.  We  note  that  it  is  easy  to  construct  smooth  generalized  “nite  element 
approximation. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OF 

18.  NUMBER 

ABSTRACT 

OF  PAGES 

40 

19a.  NAME  OF 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


first  addressed  in  [26],  and  the  term  superconvergence  was  first  used  in  [17]. 
Superconvergence  has  been  extensively  studied  ([2,  11,  18,  22,  27,  28,  33,  34,  36, 
37]  to  name  a  few)  and  there  are  more  than  1000  papers  available  on  the  subject. 
An  extensive  bibliography  (before  1998)  on  superconvergence  is  available  in 
[20],  where  as  many  references  on  3-dimensional  problems  can  be  found  in  [19]. 
Moreover,  there  have  been  several  books  written  on  superconvergence  in  the 
context  of  the  finite  element  method,  e.g .,  [1,  6,  13,  14,  23,  35,  38]. 

Typically,  superconvergence  in  the  FEM  has  been  studied  for  triangular 
meshes,  as  well  as  for  quadrilateral  meshes  with  tensor-product  elements,  but 
rarely  for  serendipity  elements.  Moreover,  the  mesh  is  required  to  have  some 
local  regularity,  e.g.,  elements  are  essentially  translation  invariant.  Also,  most 
of  these  studies  are  confined  within  the  interior  of  the  underlying  domain,  and 
only  a  few  address  the  issue  of  superconvergence  up  to  the  boundary. 

Later,  a  systematic  approach  was  introduced  in  the  analysis  of  superconver¬ 
gence  in  [9,  10],  in  the  context  of  the  finite  element  method.  This  analysis  allows 
more  general  meshes,  where  the  elements  could  be  grouped  into  translation  in¬ 
variant  “cells”  (in  contrast  to  elements  being  translation  invariant).  The  cells 
could  contain  arbitrary  number  of  elements  of  different  types.  It  was  shown  in 
these  studies  that  the  existence  of  natural  superconvergence  points  was  equiv¬ 
alent  to  the  existence  of  roots  of  a  system  of  polynomial  equations.  Moreover, 
the  superconvergence  points  are  obtained  from  these  roots,  which  (the  roots) 
are  computed  numerically.  In  special  situations,  the  system  of  equations  can  be 
written  explicitly  and  roots  can  be  computed  analytically,  as  shown  in  [36,  37]. 
In  the  context  of  finite  element  approximations  of  solutions  of  the  Poisson’s 
equation  and  the  Laplace  equation,  superconvergence  was  studied  in  [10]  for 
four  different  types  of  triangular  meshes,  as  shown  in  Figure  1.1,  as  well  as  for 
square  mesh  with  tensor-product,  intermediate,  and  serendipity  elements. 


Figure  1.1:  (a)  Regular  pattern;  (b)  Chevron  pattern;  (c)  Union  Jack 
pattern;  (d)  Cris-cross  pattern. 

We  state  some  of  these  results  (from  [9,  10];  see  also  page  354  of  [6])  in 
the  context  of  approximation  of  the  solution  of  the  Poisson  equation.  It  was 
shown  that  for  triangular  meshes,  there  are  no  natural  superconvergence  points 
for  a  mesh  (a)  with  regular  pattern  when  p  >  2  and  even;  (b)  with  Chevron 
pattern  when  p  is  even;  (c)  with  Union  Jack  pattern  for  any  p.  Also,  for  a 
square  mesh  with  serendipity  elements,  it  was  shown  that  (a)  there  are  4  natural 
superconvergence  points  and  a  superconvergence  line  for  p  =  3;  (b)  there  are 
no  natural  superconvergence  points  when  p  >  4  and  is  even,  where  as,  there 
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are  3  such  points  when  p  >  5  and  is  odd.  The  coordinates  of  all  these  points 
can  be  found  in  [10].  These  results  illustrate  the  intrinsic  complexity  of  the 
superconvergence  phenomenon,  and  they  indicate  that  the  approach  is  quite 
general  to  analyze  this  complexity.  In  [8],  the  study  of  superconvergence  was 
extended  to  the  mesh  cells  near  the  boundary  of  the  domain.  This  approach 
was  also  used  in  [7,  8]  to  study  the  effectiveness  of  various  a-posteriori  error 
estimators. 

In  this  paper,  we  will  address  the  problem  of  superconvergence  in  the  context 
of  Generalized  Finite  Element  Method  (GFEM).  This  method  was  introduced 
in  [4]  and  later  developed  and  elaborated  in  [5,  24].  It  is  a  Galerkin  method 
that  uses  a  mesh  only  minimally,  and  allows  the  use  of  non-polynomial  shape 
functions.  We  will  follow  the  approach  of  [10]  in  the  analysis  of  superconvergence 
presented  in  this  paper.  We  will  address  the  superconvergence  in  GFEM  only 
in  the  interior  of  the  underlying  domain. 

The  main  results  of  this  paper  are  Theorems  4.1  and  4.2  given  in  Section 
4.  Theorem  4.1  shows  that  the  superconvergence  points  can  be  obtained  by 
finding  the  zeros  of  a  system  of  equations  that  does  not  depend  on  the  exact 
solution  of  the  boundary  value  problem.  Theorem  4.2  addresses  the  supercon¬ 
vergence  points  for  the  second  derivatives  of  the  generalized  finite  element  error. 
We  mention  that  GFEM  allows  smooth  approximation,  in  particular  a  C 2  ap¬ 
proximation,  which  in  turn  allows  us  to  address  the  superconvergence  of  second 
derivatives  of  the  error. 

We  briefly  describe  the  organization  of  this  paper.  In  Section  2,  we  describe 
the  GFEM  and  review  the  main  approximation  result.  In  Section  3,  we  discuss 
the  so  called  interior  estimates,  which  is  crucial  for  the  superconvergence  analy¬ 
sis.  In  Section  4,  we  present  the  main  results  of  this  paper,  namely,  Theorems 
4.1  and  4.2.  We  present  an  example  in  Section  5  that  illuminates  the  results 
obtained  in  Section  4. 


2  Generalized  Finite  Element  Methods 

In  this  section,  we  briefly  describe  the  GFEM  in  the  context  of  the  approxima¬ 
tion  of  the  solution  of  a  linear  Neumann  boundary  value  problem. 

Let  12  C  R2  be  a  domain  with  piecewise  smooth  boundary  912.  We  consider 
the  Neumann  problem 

A u  =  /,  on  0, 

du  C2-1) 

—  =  g,  on  dll, 


where 


f  dx  +  /  g  ds  =  0. 


(2.2) 


ion 


We  now  give  the  standard  variational  formulation  of  the  above  problem.  Let 


B{u,v)  = 


In  l 


du  dv 
dx  dx 


du  dv 
dy  dy\ 


dx  dy 
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and 


F(v)  =  /  fvdx+  /  gvds, 

Jn  Jan 


an 

where  we  assume  /  €  L2(Q)  and  g  £  L2(dQ).  We  will  often  use  the  notation 

du  dv  du  dv 
dx  dx  dy  dy 


Bm(u,v )  =  / 

J  M 


dx  dy, 


(2.3) 


where  M  C  fi.  The  weak  formulation  of  (2.1)  reads, 

{Find  u  £  H1(  fl)  satisfying 

B{u,v)  =  F(v)  for  all  v  €  H1( fl). 

The  above  problem  is  uniquely  solvable  up  to  a  constant;  we  assume 


(2.4) 


u  dx  =  0, 


which  ensures  a  unique  solution  of  (2.4). 

The  GFEM  to  approximate  the  solution  of  (2.4)  is  a  Galerkin  method  where 
the  construction  of  trial  and  test  spaces  depend  on  a  (i)  partition  of  unity  (PU), 
and  (ii)  local  approximating  spaces. 

(i)  For  0  <  h  <  1,  a  parameter,  let  be  convex  sub-domains  of  fi 

with  N  =  N(h)  such  that  dff  =  diam(wjl)  <  2  h  for  j  =  1,2 , . . .  ,N.  We  assume 
that  for  each  value  of  h, 

N(h) 

U  <4  =  fit  (2-5) 

j=i 

and  that  any  x  £  fi  belongs  to  at  most  k  of  the  sets  where  k  is  independent 
of  h.  The  sub-domains  wj1  are  called  patches.  Clearly,  {w^l}^=1  is  an  open  cover 
of  fl.  Let  {(f>j  be  a  family  of  C 2  functions  defined  on  Q  satisfying 


&j(x,y)  =  0, 

N(h) 

X]  4>){x,y)  =  1, 
1=1 

max  \<pj(x,  y)\  <  C\, 
(x,y)en 


max  \Da(/)^(x,y)\  < 
(x,y)en 


C2 

\df\\<*\’ 


for  (x,y)  £  1  <  j  <  N(h),  (2.6) 

for  {x,y)  £  fl,  (2.7) 

for  1  <  j  <  N(h),  and  (2.8) 

for  |a|  <  2  and  1  <  j  <  N(h)  (2.9) 


where  a  is  a  multi-index  and  constants  C\ ,  C2  are  independent  of  h.  It  is  clear 
from  (2.7)  that  {4>^}jLt  form  a  partition  of  unity. 

(ii)  To  each  patch  uih ,  we  associate  an  m7-dimensional  space  Vj1  of  functions, 
defined  on  ccj1,  given  by 


Vn 

] 


^  &£•$,  b%  £  r,  $  e  H\u;^  n  C(tf) 


(2.10) 
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and  we  assume  that  V j'  contains  constant  functions.  The  space  Vj1  is  called  a 
local  approximating  space.  In  the  rest  of  this  paper,  we  will  suppress  the  h  in 
LOj,  dj,  N(h),(j>j,  V?1,  and  and  refer  to  them  as  L0j,dj,4>j,N,Vj,^j,  and 
£ji  respectively  for  notational  clarity  and  convenience,  with  the  understanding 
that  they  depend  on  h. 

The  trial  and  test  spaces  in  GFEM  is  given  by 


SGFEM  =  jy,  =  ^  where  £,•  e  Vj  | 

=  span  {//,,;  =  <f>j£ji\  1  <  i  <  rrij  and  1  <  j  <  N}  .  (2.11) 

The  functions  {pji}  are  the  shape  functions  of  SGFEM .  Finally,  the  GFEM  to 
approximate  the  solution  of  (2.4)  is  given  by 

(  Find  ugfem  G  Sgfem  satisfying 

<  fn  ugfem  =  0,  (2-12) 

[  B(ugfem,v)  =  F(v)  for  all  v  G  SGFEM . 

This  problem  has  a  unique  solution  and  is  equivalent  to  a  system  of  linear 
algebraic  equations.  Specifically,  if  we  write 

N  m.j 

ugfem  =  EE  CjiVjii 

3=1  i= 1 


then  (2.12)  yields  the  linear  system 

(  N  rrij 


iy  ■■■'J  « 

E  E  cii  /  rb*  =  ^ 

•_  1  1  J  UJn 


j  =  1  <=  1 
N  rrij 

EE5^  j  Vji)Cji  —  F{vik)i  ^  ^  —  Ulij  1  G  f  IV. 

1=1  i=l 


(2.13) 


We  note  that  the  shape  functions  {r/ji}  could  be  linearly  dependent,  and  thus 
the  dimension  of  the  null  space  of  the  matrix  in  (2.13)  could  be  greater  than 
zero.  In  this  case,  the  system  (2.13)  does  not  have  a  unique  solution.  However, 
ugfem  is  unique,  i.e.,  if  {c^}  and  {c^f  }  are  two  solutions  of  (2.13),  then 


S'1  J  L  3 

N  m, 


N  rrij 

ji 

j=  1  2—1  j=  1  2=1 


ugfem  —  EE^  =  EE^- 


For  examples  of  linearly  dependent  shape  functions  in  GFEM,  see  [3] .  In  this 
paper,  we  will  use  a  particular  SGFEM ,  which  will  have  linearly  independent 
shape  functions. 
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We  now  present  two  results  on  approximation  properties  of  SGFEM ,  which 
in  turn  give  an  error  estimate  for  ||it  —  ugfem\\h1(q.)-  Suppose  u  £  H1( Q)  can 
be  accurately  approximated  on  LOj  by  a  function  £  Vf,  specifically,  suppose 


\u-tt 


<  4(j) 


and 


Define 


l«-£l  h>m<4U)- 

N 

Cu  \  ^  ±  cu  r~  qGF EM 

Z  =  E  <!>&  G  s 

j=i 


Then  we  have  the  following  two  results  ([5,  24]). 


Theorem  2.1  Suppose  u  £  H1{  12).  Then 


n«-r 


|l2(0) 


<  C\K1/2 


N 

E<? 


1/2 


(j) 


and, 


N 


1/2 


It* - fV(n)  <  (2«)1/2  E^jjr  +  ^E 4(3) 

where  C\,  Ci  are  as  in  (2.8),  (2.9),  respectively,  and  dj  =  diam(u>j).  □ 


Remark  2.1  We  note  that  the  above  theorem  is  true  even  when  the  patches 
LOj  are  non-convex. 


Theorem  2.2  Suppose  u  £  Hl(£i).  Suppose  the  patches  LOj  satisfy  the  following 
assumption: 

For  all  1  <  j  <  N,  there  exists  C3  >  0,  independent  of  j,  such  that 
\Ml2(u }j)  <  C3dj  \v\Hi(u jj),  for  all  v  G  H1( uj)  satisfying  /  vdx  =  0,  (2.14) 

Jujj 

where  dj  =  diam(u>j). 

Then,  there  exists  £  V)  such  that 


N 

E  °jZ'j  e  SGFEM 
1=1 


satisfies 


llu  -  l“IU2(n)  <  C5 


1/2 
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and 


1/2 


N 


\u  ~  £“|ff1(^)  <  C6  (  ^2  el (j) 
<j=1 


where  C$,Cq  depend  on  C\,  C2,  C^.D 


Remark  2.2  It  is  shown  in  [3]  that  (2.14)  holds  when  the  patches  u)j  are  convex. 
The  precise  dependence  of  C5 ,  C%  on  C'i ,  C  2 ,  C3 ,  and  the  dependence  of  the 
Poincare  constant  C3  on  the  geometric  data  of  tOj  is  also  given  in  [3] .  In  the  rest 
of  the  paper,  we  will  not  differentiate  between  various  constants,  and  instead 
will  use  a  generic  constant  C. 


Theorem  2.2  gives  an  error  estimate  for  the  GFEM.  Suppose  the  hypothesis 
(2.14)  is  satisfied  and  suppose  u  is  the  solution  of  (2.4).  Then  from  Theorem 
2.2  we  have 


N 


1/2 


I|w  —  UgfEmWh !(Q)  <  ||u  —  <  C  £  4ti) 


U=1 


It  will  be  useful  to  state  this  estimate  in  the  form 


II  u 


UGFEm\\h1(Q.)  <  C 


1/2 


\w~i 


(2.15) 


To  obtain  the  main  result  of  this  paper,  we  will  impose  additional  restrictions 
on  the  patches  {w^},  the  partition  of  unity  {<(/,■},  and  the  local  approximation 
spaces  Vj.  We  list  them  as  three  assumptions. 

Assumption  A  In  addition  to  (2.5),  we  assume  that 


UjCUj,  1  <j<N, 

where  to*  is  a  ball  of  diameter  d*,  and  there  exists  0  <  er  <  1,  independent  of 
the  parameter  h,  such  that 


d*>adj,  1  <j<N.  (2.16) 

Assumption  B  In  addition  to  (2.6)-(2.9),  we  assume  that 
<t>j(x,y)  =  1,  for  (x,  y)  G  to* . 

Since  {4>j}  is  a  partition  of  unity,  it  is  clear  that  to*  fl  w*  =  0  for  j  ^  *,  and 
4>j(x,  y)  =  0  for  ( x ,  y)  €  to*  when  i  ^  j. 

Assumption  C  We  consider  V7  =  1  <  j  <  N,  where  Vk{tOj)  is  the 

space  of  polynomials  of  degree  k  on  tOj.  We  assume  that  for  1  <  p  <  00, 

(a)  IICIIw^)  <  C||^llwps(wp;  for  all  £  G  V)  and  0  <  s  <  k,  (2.17) 
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where  C  depends  on  k,  but  not  on  j  or  h: 

(b)  1 1  £  1 1  Wp  (ujj )  <  CdJ{t-s)U\\w^),  for  all  £  G  Vj  and  0  <  s  <  t  <  k.  (2.18) 

where  C  depends  on  s  and  t,  but  is  independent  of  j  and  h. 

Furthermore,  when  ujj  0  dfl  7^  0,  we  may  allow  to*  to  be  a  portion  of  a  ball 
bj  (jf  O  with  center  inside  u>j ,  such  that  lo*  =  bj  C l  f2  satisfies  uj*  C  u)j  and 

0<B1<yi<B2<l,  (2.19) 

I  yj  I 


where  the  constants  B\,B2  are  independent  of  j  and  h ,  and  |w*|,  \bj\  are  the 
areas  of  the  sets  u*,  bj  respectively. 

We  note  that  (2.17)  holds  if  the  patches  LUj  satisfy  some  reasonable  assump¬ 
tions.  For  example,  let  Bj  be  the  smallest  ball  containing  ojj  with  the  same 
center  as  ui*.  If 


x  <  diam(Bj)  <  ^ 
—  diam(oj*)  ~ 


(2.20) 


then  one  can  show  that  (2.17)  is  satisfied. 


We  will  often  denote  SGFEM  by 

Sh  =  Sh’k,  (2.21) 

where  k  indicates  degree  of  the  polynomials  used  in  the  Vj  s.  Let  {ijiY^d-y  be  a 
basis  of  Vj  for  j  =  1, 2, . . . ,  N.  Then 


Vji  =  1  <  *  <  m.j,  1  <  j  <  N 

are  the  shape  functions  of  Sh.  In  the  following  proposition,  we  give  an  easy 
proof  that  {rjji}  is  a  linearly  independent  set. 


Proposition  2.1  The  set  of  shape  functions  {i)ji,  1  <  i  <  m.j,  1  <  j  <  N}  is 
linearly  independent. 


Proof:  Suppose  the  set  is  linearly  dependent  and  there  are  constants  c?!;,  1  < 
i  <  m.j ,  1  <  j  <  N,  not  all  zero  such  that 

N  rrij 

EE  CjiT]ji(x,y )  =  0,  for  all  ( x,y )  €  O.  (2.22) 

i= 1  *= 1 

Without  loss  of  generality,  suppose  the  constant  c,0.j0  7^  0  for  some  1  <  io  < 
mj,  1  <  jo  <  N.  From  Assumption  B,  we  have  (j>j0(x,y)  =  1  for  ( x,y )  £  lo*o, 
and  < pj(x,y )  =  0  for  (x,y)  G  uj*o,  j  7^  j0-  Therefore  using  (2.22),  we  get 

N  mj  mjQ 

EEr"  'lji(x.y)  =  ^ ~2cj0ii(t>j0(x,y)^j0ti(x,y ) 

j—1  i—1  i—1 

m30 

=  E cio,£ioAx’ v)  =  for  a11  (®) y )  e  oj*0. 

i= 1 
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But  {£j0j}i=?°  is  a  basis  of  Vj0,  and  thus  we  conclude  from  above  that  Cj0>l  =  0 
for  1  <  i  <  m.j0.  In  particular,  we  have  Cj0^0  =  0,  which  is  a  contradiction. 
Hence  {rjji,  1  <  *  <  nij,  1  <  j  <  N}  is  linearly  independent.  □ 

As  mentioned  before,  the  set  of  shape  functions  {fjrl}  for  the  GFEM  may 
be  linearly  dependent.  In  contrast,  the  set  of  shape  functions  for  the  GFEM 
obtained  using  the  patches  ojj  and  the  partition  of  unity  functions  <f>j  satisfying 
assumptions  A  and  B  is  linearly  independent.  Thus  the  linear  system  (2.13)  has 
a  unique  solution. 

We  now  present  a  few  examples  of  patches  {toj}  and  partition  of  unity  func¬ 
tions,  defined  relative  to  these  patches,  satisfying  assumptions  A  and  B. 


Example  1:  Let  a  €  R  be  such  that  0  <  er  <  1.  Consider  r  G  K.  such  that 

For  r  fixed,  let  s(a;)  be  a  smooth  function  on  the  interval  [r,  1  —  r]  satisfying 


s(r)  =  1,  s(l  —  r)  =  0, 

s^\r)  =  s(t)(l  —  r)  =  0,  for  t  =  1,  2, . . . ,  l 


We  now  define  a  smooth  function  <p(x)  on  [—1, 1]  by 

\x\  <  r 
M  >  1  —  r 
r  <  x  <  1  —  r 
—  (1  —  r)  <  x  <  —r 

Clearly,  c/)(x )  G  Cl(— 1,1)  and  support  of  c/)(x)  is  [—(1  —  4,(1  —  r)]-  We  also 
note  that 


<t>(x )  = 


(j){x)  +  (f){x  —  1)  =  s(x)  +  1  —  s(x)  =  1,  for  r  <  x  <  1  —  r.  (2.24) 

Suppose  H  =  (0, 1)  and  we  consider  the  nodes  Xi  =  ih,  i  =  0,1,2 ,N, 
where  Nh  =  1.  For  *  =  1, 2, . . , ,  N  —  1,  we  define  patches 

w,:  =  ( Xi  -  (1  -  r)h,  Xi  +  (1  -  r)h) ,  (2.25) 

to*  =  (xi  —  rh ,  Xi  +  rh ) 

For  i  =  0,  N,  we  define 

wo  =  (0,  (1  -  r)h ) ,  u)q  =  (0,  rh) 
u>N  =  (1  -  (1  -  r)h,  1) ,  =  (1  -  rh,  1)  (2.26) 


Clearly,  =  G  and  to*  C  tOi  for  i  =  0, 1, 2, . . . ,  N.  Also  using  (2.23),  we 

can  easily  show  that 

-j-  =  ——  >  a,  for  i  =  0, 1, . . . ,  N, 

di  1  -  r 
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where  eft  =  diam(wi)  and  d*  =  diamftu*).  Thus  the  patches  {cu,}  satisfy  As¬ 
sumption  A. 

Now  for  each  i,  0  <  i  <  N,  define  functions  <f>i(x)  on  12  by 


<pi(x) 


41  (—f f1)  ,  Z  G  [a:*  -  ft,  Xi  +  ft] 
0,  x  G  12  \[a’j  —  ft,  Xi  +  ft] 


(2.27) 


It  is  easy  to  check  that  <j>i(x)  =  0  for  ie!l  \w»,  and  <j)i(x)  =  1  for  x  G  u>*  (which 
is  Assumption  B).  Also  from  a  standard  scaling  argument,  it  is  immediate  that 


\<t>i\x)\  < 


Ct 

|diam(w,;)|t  ’ 


0  <  t  <  ft 


where  Ct  depends  only  on  max^^n  \^f\y)\- 

We  now  show  that  form  a  partition  of  unity.  Consider  the  interval 


[xi,  xi+\]  =  [xi,  Xi  +  rft]  U  [xi  +  rh ,  £i  +  (1  -  r)h\  U  [xi  +  (1  -  r)ft,  xi+\}. 

For  x  G  [a;*,  a;,  +  rft],  we  have  0  <  (x  —  Xi)/h  <  r,  and  from  the  definition  of  <j>i 
and  <f>, 

4ft  0*0  =  Mx)  =  (/>  (  ,  Xl  )  =  1- 

Similarly,  for  xG  [xi  +  (1  —  r)ft,  Xi+i], 


N 

X]  4ft  0*0  =  <t>i+l{x)  =  1. 

0=0 


For  x  G  [a^  +  r/i,  +  (1  —  r)/i],  we  have  r  <  (x  —  Xj)/h  <  1  —  r,  and  from  (2.24), 
N 

j=0 


Thus 

N 

'y^  <j)j(x)  =  1,  for  all  x  G  11. 
i— o 

We  note  that  for  a  two  dimensional  domain  =  (0, 1)  x  (0, 1),  it  is  possible 
to  construct  patches  of  the  form  Ui  x  u>j  and  the  partition  of  unity  function  of 
the  form  4>i(x)(j>j(y)  that  satisfy  Assumptions  A  and  B.  We  do  not  describe  this 
construction  in  detail  here. 

Example  2:  Let  12  be  a  domain  in  R2.  For  0  <  ft  <  1,  we  consider  the  points 
{xj  =  (Xj,yj)}jL1  such  that  Xj  €  12.  We  will  suppress  ft  in  Xj,  and  instead 
will  denote  these  points  by  xr  We  assume  that  the  points  are  distributed  in  a 
way  such  that  the  following  hold: 


<t>i{x)  +  Oi  i.i(a-) 
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(a)  For  each  j,  there  is  a  disc  O*  of  radius  r*  =  O(h)  and  centered  at  x:)  such 
O*  (~l  O*  =  0  for  i  ^  j.  Let  to*  be  the  disc  centered  at  Xj  of  radius  r*  / 2. 
If  w*  (jL  0,  we  consider  to**  instead  of  to*,  where  to**  =  to*  fl  O  satisfies 
(2.18)  with  bj  =  to*  and  to*  =  to**.  We  redefine  to*  as  to**. 

(b)  For  each  j,  there  is  a  convex  open  set  (patch)  tOj  with  diam(wj)  =  dj  = 
O(h)  such  that  O*  C  to  j  and  there  exists  0  <  a  <  1  such  that 

r*  >  adj ,  1  <  j  <  N. 

Moreover,  {tOj}j r=1  form  an  open  cover  of  fl  satisfying  (2.4).  It  is  easy  to 
check  that  Assumption  A  is  satisfied. 


We  will  now  construct  partition  of  unity  functions  subordinate  to  the 

covering  {tOj},  satisfying  Assumption  B.  For  each  j,  we  first  consider  a  smooth 
non-negative  function  0  <  ipj(x,  y)  <  1  on  f 1  such  that 


0,  (a;,  y)  G  to* 

1,  (x,y)  G  Q\0* 


(2.28) 


(ii)  max  \Daipj\  <  C/h)a\,  for  |a|  <  2.  (2.29) 

Men 

The  function  tpj  could  be  a  radial  function  based  on  a  scaled  and  suitably  defined 
one-dimensional  function  in  [0,oo).  We  then  define  the  function 

N 

^(x,y)  =  Y[i/Jj(x,y). 
i 

Clearly,  0  <  ip(x,y)  <  1  and  using  (2.28),  (2.29),  we  can  show  that 


f  °. 

(x,y)  G  to*,  1  <  j  <  N 

i>(%,y)  =  <  !. 

(x,  y)  G  0\  U f=10* 

(2.30) 

1  ’•Pjix^y), 

(x,y)  G  O*,  1  <  j  <  N 

max  \Dail)\  <  C/h)a\ 
(x,y)  6f2 

,  for  a  <  2 

(2.31) 

We  next  consider  smooth  non-negative  functions  fj  (x,  y)  with  compact  sup¬ 
port  in  the  patch  tOj,  satisfying  max(IiV)£w  \Dafj\  <  C/h)a '  for  |a|  <  2.  The 
functions  fj(x,y)  could  be  constructed  as  radial  functions  with  circular  support. 
Also,  the  construction  of  functions  fj(x,y)  with  polygonal  support  have  been 
discussed  in  [16,  15].  We  further  assume  that  there  exists  7  >  0  such  that 

N 

'^2fj(x,y)>'j>  0,  for  (x, y)  G  fi,  (2.32) 

3= 1 

and  define 

y)  =  fj(x, y)ip(x ,  y)  +  [1  -  ipj(x, y)],  1  <j<  n. 

We  now  state  some  relevant  properties  of  <j)j(x,  y). 
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(i)  The  supp  <j>j{x,y)  C  LOj. 

This  is  clear  from  the  fact  that  supp fj(x,y)  C  uoj  and  il>j(x,y)  =  1  for 
(. x,y )  e  £2\0*  (see  (2.28). 

(ii)  <i>j(x,  y)  =  Sij  for  (x,y)  Gw*  ). 

For  (a:,  y)  Gw*,  we  have  from  (2.28)  and  (2.30),  that  ipj(x,  y)  =  y)  =  0 
and  thus  <j>j(x,  y)  =  1.  Similarly,  we  show  that  4>j(x,  y)  =  0  for  (x,  y)  G  w*, 
i  ±  3- 

(iii)  max(Xi„)eu.  \Dacj)j\  <  C/hH,  for  |a|  <  2. 

This  is  obtained  using  (2.29),  (2.31),  and  the  fact  that  max^^g^.  \Dafj\  < 
C/h l“l  for  \a\  <  2. 

(iv)  There  is  7  >  0  such  that  <j>j(x,y)  >  7,  for  (x,y)  G  £2. 

To  obtain  this  result,  we  first  note  from  (ii)  that  for  1  <  i  <  N  and 
for  (x,y)  G  w*,  we  have  <j>i(x,y)  =  1  and  <j>j(x,y)  =  0  for  j  ^  i.  Thus 

Y^=ihix^v)  =  Mx,y)  = 1- 

Also  for  (x,  y)  G  f2\U^1  O*,  we  have  from  (2.28)  and  (2.30)  that  ip(x,  y)  = 
1,  ipj(x,y)  =  1  for  all  j,  and  thus  from  (2.32)  we  get  ]T7=1  4>ji.xiy)  >  7- 

Moreover,  for  (. x,y )  G  0*\w*,  1  <  i  <  N,  we  have  ipj(x,y)  =  1  for  j  ^  i 
and  consequently,  using  (2.32)  and  assuming  7  <  1,  we  get 

N 

^2^j{x,y)  >  ii>i(x,y)  +  l-il)i(x,y) 

3- 1 

>  1  -  1(7- i)IV’i(a;,y)  >  1 +  7- 1  =  7- 

If  7  >  1,  it  is  easy  to  show  that  4>j(x^y)  >  1  f°r  (xi  y)  G  0*\w*, 

1  <  i  <  N. 

Thus  the  result  is  true  with  7  =  min  (1,7). 


Finally,  we  use  the  technique  of  Shepard  ([21,  31])  to  define 


<t>j{x,y) 


4>jix,y) 

T,^LiMx,y) 


Using  the  properties  (i)-(iv)  of  4)j(x,  y),  given  above,  it  is  easy  to  check  that 
{(j)j}jL  1  is  a  partition  of  unity  subordinate  to  the  cover  satisfying  (2.6)- 

(2.9)  and  Assumption  B. 

The  points  {xj}!'[=1  considered  in  this  example  have  to  be  distributed  such 
that  (a)  and  (b),  mentioned  above,  are  satisfied.  For  example,  if  the  particles 
Xj  are  vertices  of  a  quasi-uniform  triangulation  of  £2,  it  is  possible  to  construct 
w*,  O* ,  and  uij  satisfying  (a)  and  (b). 
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3  Interior  Estimate 


Interior  estimates  play  a  crucial  role  in  the  study  of  superconvergence  in  a 
Galerkin  method.  In  a  series  of  papers  ([25,  29,  30]),  Nitsche,  Schatz,  and 
Wahlbin  developed  a  machinery  to  establish  interior  estimates  in  the  context  of 
finite  element  method.  This  theory,  which  is  based  on  certain  axioms  on  the 
finite  dimensional  approximating  subspace,  can  also  be  used  in  the  context  of 
GFEM.  In  this  section,  we  will  show  that  the  finite  dimensional  space  used  in 
GFEM,  i.e.,  Sh  =  SGFEM ,  satisfies  the  axioms  given  in  [30]. 

We  first  state  the  interior  estimate  that  we  will  use  in  the  next  section.  Let 
fin  CC  f Id  CC  fi  be  domains,  where  D  =  dist(flo,  dCln)-  We  also  assume  that 
all  the  patches  u>i  s  in  a  neighborhood  of  are  quasi-uniform,  i.e.,  0  <  /3  < 
di/h.  We  also  assume  that  D  >  cq h  for  co  large  enough.  Let  Uh  €  Sh(flD)  = 
SGFEM (flD)  be  such  that 

B{u  —  Uh,v)  =  0,  for  all  v  €  Sh(flD)-  (3.1) 

Here  Sh(Q]j)  denotes  the  restrictions  of  functions  in  Sh(fl)  to  f Ip,  and  Sh(Q£,) 
denotes  the  set  of  functions  in  Sh(£li))  with  compact  support  in  the  interior  of 
Qd-  We  now  state  Theorem  1.2  from  [30],  which  will  be  used  later  in  this  paper. 

Theorem  3.1  (Theorem  1.2  of  [30])  There  exists  a  constant  C,  depending 
only  on  the  constants  in  the  Axioms  A1-A5  (given  below)  over  fljj,  such  that  if 
e  =  u  —  Uh  satisfies  (3.1),  then 

\e\ w^(q.0)  +  D  1||e||Loo(o0) 

<  C  min  (|u  -  xl w^irin)  +  D~x\\u  -  xIIwod)) 

X&sh  ” 

+CD-2||e||i2(nD)).D  (3.2) 

Theorem  1.2  of  [30]  is  quite  general.  The  theorem,  stated  above,  can  be  obtained 
by  using  s  =  0,  q  =  2,  and  the  fact  that  fi  C  K2  in  Theorem  1.2  of  [30]  . 

The  above  theorem  holds  provided  the  subspace  Sh(fl)  satisfies  certain  ax¬ 
ioms.  For  G  C  fi,  let  Sh(G)  be  the  restriction  of  Sh(il)  to  G,  and  let 

Sh(G)  =  {X  :  X  e  Sft  (G),  supp  x  CC  G}. 

Also,  for  A  CC  fl,  define 

7(A)  =  {jGN  :  A  fl  0} 

and 

A=  (J  Wj. 

ie  7(A) 

It  is  clear  that  A  c  A. 
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We  now  show  that  there  exists  positive  constants  Gi,  C2,  G3,  C4,  /co,  and 
ho  <  1  such  that  the  space  Sh  =  SGFEM  satisfies  the  following  axioms. 

Axiom  Al.  Approximation.  Let  G2  C  G  CC  O  with  dist(G,  dfl)  >  k$h. 
Then  for  each  v  €  Wlq(G2),  there  exists  x  €  Sh(G)  such  that  for  G\  CC  G2 
with  dist(Gi,  9G2)  >  koh, 

11^  ~  xllw'(Gi)  <  Cih1  ‘Hull w‘(g2)  (3-3) 

for  0  <  t  <  l  <  k  +  1,  1  <  q  <  00,  t  =  0, 1, 2. 

Moreover,  if  supp  v  C  Gi,  then  x  €  5lft  (G2). 

Remark  3.1  We  note  that  in  [30],  the  Wq’h(Gi)  norm  was  used  in  (3.3)  instead 
of  W*(Gi).  It  was  natural  to  use  Wq,h(G\)  norm  in  [30],  since  the  space  Sh 
considered  in  [30]  was  a  subset  of  fl  C2'h,  i.e.,  the  functions  is  Sh  were 
piecewise  G2  and  globally  functions.  We  note  that  functions  x  £  S1*  = 
gGFEM  are  (j2  functions,  and  thus  Hxllw^GGi)  =  llxllw|(Gi)-  For  a  definition 
of  Wq,h(G\)  norm,  we  refer  to  page  925  of  [30]. 

Proposition  3.1  The  subspace  Sh  =  SGFEM  satisfies  Axiom  Al. 

Proof:  Let  v  be  an  extension  of  v  to  fi,  such  that  v  =  v  in  G2  and 

ll^ll^(Q)  <  C1MI  w‘(G2)-  (3-4) 

For  existence  of  such  an  extension,  we  refer  to  [32] . 

Let  Gf  be  the  averaged  Taylor  polynomial  of  v  of  degree  l  —  1,  averaged  over 
to*.  Then  Gf  \ _  £  Vj,  and  from  Lemma  4.3.8  of  [12],  we  know  that 


<  I 


I  WUu 


Define 


Clearly,  x  €  Sh(Cl)  and 


1= 1 


(3.5) 


N 


K’  -  xl 


w*m 


<cY,ofir  Gi 


\q 

1  WUu 


1=1 
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where  G  depends  on  k  and  q.  Therefore,  using  (2.9),  (3.5),  and  (3.4),  we  get 


\v~x\q 


Hence, 


< 


wnm  - 


< 


v\ q 


N  t 

j= 1  i=0 
N  t 


j— 1  i= 0 
N 


< 


j= i 


W'M 


<  ChV-»q\v\ll{Q)  <  ChV-** 


1 wug2 


lW  Xlw^(Gi)  —  \v  Xl  W*(Q)  —  t)q\\v\\w^(G2), 

from  which  we  get 

llv  -  x\\w*(G$  —  Ch*'1  t'>9\\v\\wUG2)' 


Now  letting  %  =  x| G  €  Sh(G)  we  get  the  (3.3). 

We  now  suppose  that  supp  v  C  G\ .  We  consider  ko  large  enough  that 
G\  CC  G2.  Since  supp  v  C  Gi,  we  have  =  0  for  j  ^  7(Gi),  where  v  is 

defined  as  the  zero  extension  of  v.  Thus  supp  x  C  Gi  CC  G 2,  and  hence, 
X  €  Sh(G2),  which  proves  the  desired  result.  □ 


Axiom  A2.  Inverse  inequality.  Let  G 1  CC  G2  with  dist(Gi, 9G2)  > 
koh.  Then  for  x  €  Sh(G2), 

llxllffMGi)  <  G2/i_1||x||l2(G2)-  (3-6) 

Moreover, 

||xl|w*(Gi)  <  C2ht-s-2^-1^\\x\\wk(G2),  (3.7) 

for  0  <  t  <  s  <  2,  1  <  qi  <  q  <  00. 


Remark  3.2  We  note  that  (3.6)  is  a  special  case  of  a  more  general  inverse 
inequality  assumption  given  in  [30] .  A  careful  reading  of  the  proof  of  Theorem 
1.2  in  [30]  shows  that  we  need  only  (3.6)  to  get  Theorem  3.1  in  this  paper. 


Proposition  3.2  The  subspace  Sh  =  SGFEM  satisfies  Axiom  A2. 
Proof:  Suppose  x  C  Sh(G2)-  Then 

x\Gl  =  where  &  €  Vi. 

*S7(Gi) 


15 


Recalling  that  fa  =  0  in  f Y\u>i,  we  have 

MhRGi)  -  ^  E 

iey(Gi) 

=  C  E 

iEy(Gi) 

^  C  E  (l&lw^wdll&lli^Wi)  +  ll'AilllooCwijI^I/fRwi))  » 

i67(Gi) 


where  C  depends  only  on  k.  Therefore  from  (2.8),  (2.9),  and  (2.18)  with  p  =  2, 
we  have, 

IxIffRG i)  —  ^  2  E  ll^;l|22(wi)-  (3-8) 

*G7(Gi) 

We  now  consider  the  domain  A  such  that  G \  CC  i  CC  (?2  for  fco  sufficiently 
large.  Then  from  (2.17),  we  get 

llxlli,(.4)  >  llxlliio.)  >  E  IIxIILk) 

i£7(Gi) 

=  E  INILk)>c  E  H^llL(^)’ 

j£7(Gi)  *67(Gi) 

and  thus  from  (3.8),  we  have 

IIxIIhrgi)  <  C'^_1IIxIU2(a)  ^  Ch  1IIxIU2(g2), 

which  is  (3.6). 

We  now  prove  (3.7)  for  s  =  2  and  t  =  0.  The  other  cases  can  be  proved 
similarly.  Using  the  argument  leading  to  (3.8),  but  employing  (2.18)  with  p  =  q, 
and  the  fact  that 


lie, 


\Lq(ui 


<  C(h2)1/q~1/gi  ||£j| 


Lqi  (coi 


we  get, 


lx|L,(Cl,  <  Ch~*>  E  ll<< 
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*£7(Gi) 

<  Ch~2q(h2)1~q^qi  E  ll€illL(U4)* 


*67(Gi) 


Therefore, 


lxlw„2(G!)  <  Ch  ( h  ) 


-2('^2U/9-1/91 
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and  since  q±  <  q ,  we  have 


1/91 


ixi^(Gl)  <  ch-ih^-'M  m 


91 

Lqi(ui) 


U£7(Gi 


Again  as  before,  but  using  (2.17)  with  p  =  qi,  we  can  show  that 


I  Lgi(G2)-C  22  11^*11  Lqi(u>i)’ 

iej(Gi) 


and  thus  from  (3.9)  we  get 

lx|w|(Gl)  <  Ch-2~2Wqi-1/q)\\x\\Lgi(G2) 

Similarly  we  can  show  that 

lxlw|(Gl)  ^  Ch~J~2(1/qi~1/q)  ||x||l9i  (G2),  for  j  =  0,1, 

and  therefore  using  (3.10)  we  get 

||x||w|(Gl)  <  Ch~2~2<'1/gt"1/q)\\x\\Lqi(G2), 


(3.9) 


(3.10) 


which  is  the  desired  result.  □ 

Axiom  A3.  Superapproximation.  Let  G\  CC  G 2  CC  G3  with  dist(Gi,  9G2)  > 
koh,  dist(G2, dG3)  >  koh  and  let  p  £  G°°(Gi).  Then  for  each  x  €  Sh(G3),  there 
exists  an  77  £  Sh(G3 )  such  that  for  some  t  >  0, 


\\PX~v\\h*(g3)  <  C3h\\p\\weoiG1)\\x\\H^G2),  s  =  0,l,  (3.11) 

and 

IMIl„(G3)  <  C,||xl|i,(G3)>  l<9<oo.  (3.12) 

Furthermore,  let  G_ 2  CC  G_  1  CC  Go  CC  Gi  with  dist(G_2,  <9G_i)  >  fco/t, 
dist(G_i,  9Go)  >  koh.  and  dist(Go,9Gi)  >  koh.  Then,  if  p  =  1  on  Go,  we  have 
rj  =  x  on  G_i  and 


II PX  -  vWh^Gs)  <  G3/i||x||ff3(G3\G_2).  (3.13) 

We  first  prove  the  following  lemma. 

Lemma  3.1  Let  p  he  a  smooth  function  on  Ui  and  £  V).  Then  there  exists 
£  Vi  such  that 

l|P&  -  <  C'^l|p||w^(wi)||6l|fl's(w«)j  s  =  0,  1.  (3.14) 
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Proof:  We  prove  (3.14)  for  s  =  1.  The  proof  for  s  =  0  is  similar. 

Recalling  that  V)  contains  constants,  we  first  decompose  fi  as 

fi  =  fil  +  fi2, 

where  fa,  ^2  G  V)  and  £,;l  is  a  constant  that  is  orthogonal  to  fa  in  -o(wj)  inner 
product,  ie.,  <  ^1,^2  >//1(Wi)=  0-  Clearly,  fa  dx  =  0  and  from  (2.14),  we 
have 

||&2||z,2(u>j)  <  (3.15) 

We  will  now  construct  fi  =  +  fa  G  V):  such  that  (3.14)  holds  with  s  =  1. 

Let  L  G  Vl(u)i)  C  V)  such  that 

Up  -  M\ w^(ui)  <  (3.16) 

(L  could  be  the  linear  Taylor  polynomial  of  p  centered  at  the  center  of  u>*).  We 
choose  fa  =  fi\L.  Clearly,  fa  G  Vi,  and  using  (3.16),  we  have 

\\pfn  -  fnWmiui)  =  ||£a(p  -  -^)llffi(wi) 

<  C\\p  -  L\\wZa(wi)\\€il\\L2(wi) 

<  Ch\p\W2a(ni)\\fii\\m(ui)-  (3-17) 

We  next  write  p(x)  =  p+p*(x),  where  p  is  a  constant  (we  may  take  p  =  p{xf), 
where  Xi  is  the  center  of  lo*  ) ,  and 

\\p*\\LM<Ch\\Vp\\LM.  (3.18) 

We  now  choose  fa  =  pfa-  Clearly,  fa  G  V)  and  using  (3.15)  and  (3.18),  we 
have 

\\pfi2-fa\\m(uH)  =  \\p*fi2\\m(u>i) 

<  Clip*  lUoo(wi)  11^*2  ||ffi(Wj)  +  C'||Vp|U00(&)i)||^2||i2(W() 

<  Ch\\Wp\\Loo(wi)\\fa\\H,(m).  (3.19) 

Finally,  using  (3.17),  (3.19),  and  the  fact  that  <  fa,  fa  >h1( Wi)  =  0,  we  have 
\\pfi  ~  fi\\m(u>i)  <  ||p6l  -  &l||.ff1(u;4)  +  1 1 P^i2  -  fi2\\m(wi) 

<  C/l|p|y^(w.)||£,:i||.H-1(u;i)  +  Ch\\Vp\\Loo(ui)\\fi2\\HHui) 

<  Ch\\p\\w^u.)\\fi\\Hx(u.), 
which  is  the  desired  result.  □ 

Proposition  3.3  The  subspace  Sh  =  SGFEM  satisfies  Axiom  AS. 

Proof:  Suppose  \  €  Sh(G3).  Since  p  G  C°°(G  1),  the  function  p\  has  the 
form 

PX  =  ^2  P^*’  where  ft  G  V. 

i67(Gi) 
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Let  Pi{p£i)  £  Vi  be  the  ^-projection  of  pf;  onto  V,  for  i  £  7(Gi)-  Then  by 
Lemma  3.1,  we  have 


\\pii~  Pi(pii)\\m{ui)  <  \\p£,i-£i\\m(ui)  <  ChWpWw^iw^MiWH^uji),  *  €  l(Gi)- 

(3.20) 

Also  since  P,(p£,)  is  the  H 1  projection  onto  V.t,  and  V)  contains  constants,  we 
have  f  [p^i  —  Pi{pii)]  dx  =  0,  we  have  from  (2.14), 

II pti  ~  PiipmL^i)  <  Ch\pii  -  PMi) \mM.  (3.21) 

We  now  define 

v=  E 

i67(Gi) 

Clearly,  77  £  Sh(G3)  and  p  =  0  in  fl\Gi-  We  choose  k0  large  enough  such  that 
G\  CC  Gi.  Thus  77  £  Sh(G3).  Now 

IIpx-^IIhi(g3)  =  IIpx  —  7?II^i(gi) 

<  c  e  UMi-PMiK^) 

iS7(Gi) 

=  c  e  WMi-Piipmm^) 

iS7(Gi) 

<  c1  E  (ii^iiL(Wi)i^-p*(^)i^(^) 

iS7(Gi) 

+  l|V0i||L(a,i)H^  -  PMJWU^))  (3-22) 

Thus  using  (2.9),  (3.21),  and  (3.20)  in  the  above  inequality  and  then  using 
(2.17),  we  get 

IIPX-??IIhi(G3)  <  c  E  l?&  ~  Pi(P&)\2HVui) 

iS7(Gi) 

<  C'/i2||p||^2)(Gi)  E  ll^llfl-qwj) 

iS7(Gi) 

<  C'/i2||p||rv23(Gi)  E  H&llffbO'  (3.23) 

i67(Gi) 

Since  Gi  C  G2,  we  get 

IIxIIhi(g2)  —  E  IIxIIhi(o>*)  =  E  ll&llWr)’ 

*67(Gi)  *67(Gi) 

and  combining  it  with  (3.23)  we  have 

II PX  -  v\\m(G3)  <  C,^l|p|l^(G1)llxl|ffi(G2) 
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which  is  (3.11)  for  s  =  1.  The  case  for  s  =  0  can  be  proved  similarly. 

We  now  prove  (3.12).  We  first  note  from  (3.20)  and  (2.18)  that 

||-Pi(p&)IU2(wi)  <  \\p£,i\\L2(wi)  +  ChWpWw^^WiiWm^i) 

<  IMUoo(wi)ll&IU2(wi)  +  c'IICillL2(wi) 

<  c\mL2(w& 

Thus  using  (2.17)  and  the  above,  we  get 

IMIL(G3)  <  ^  E  \\<t>ipMi)\\l^) 

*67(Gi) 

<  c  ^  \\PMi)\\UUi) 

iej(Gi) 

<  c  E  n^iiL(^) 

<£7(Gi) 

<  C  E  II&IILk)  (3-24) 

<£7(Gi) 


Now, 

llxlli2(G3)  >  E  llxlll2K)=  E  H&IILkw 

»£7(Gi)  »£7(Gi) 

and  therefore  from  (3.24),  we  have 

IMU2(g3)  <  C'IIxIIl2(g3)- 

Finally,  using  the  fact  that  Sh(G3 )  is  finite  dimensional,  we  get 
\\V\ |l9(G3)  <  C||xll l9(G3) i  for  1  <  q  <  oo. 

We  now  prove  (3.13).  We  first  assume  that  dist(G_i, dG0)  >  k0h  for  a 
suitable  fco  such  that  G_i  C  Go-  Since  p  =  1  on  Go,  we  have  p  =  1  on  07 
for  i  £  7(G_i).  Also  from  the  definition  of  if1  projection,  we  have  Pi(p£,i)  = 
Pi(£i)  =  £*  for  *  e  7(G_i).  Therefore,  for  a;  €  G_i, 

f?(z)  =  E  <t>i(x)pi(p£i)(x)  =  E  fcW&ix)  =  X(x)- 

*£7(Gi)  *£7(Gi) 


Thus  rj  =  x  on  G_i. 

Now  using  the  argument  leading  to  (3.22)  and  using  (3.21)  and  (2.9)  we 
have, 

IIPX  77llli'1(G3)  —  ^  E  \P€i  ~  Pi(p£i)\'H1(ui)‘  (3.25) 

i£7(Gi) 
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We  note  that  7(G_i)  C  l(Gi).  Also  recall  that  p  =  1  and  Pi{p^i)  =  &  on  ojj 
for  i  G  7(G_i).  Thus  from  (3.25),  (3.20),  and  (2.17),  we  get 


< 

*  E 

1  Pfi  -  Pi(pfi 

)  Iff(wi) 

*67(Gi)\7(G_i 

) 

< 

IIP  IvV^(Gi) 

E 

11^11^(0,0 

*67(Gi)\7(G_i) 

< 

\\P  IvV^(Gi) 

E 

ie7(Gi)\7(G_i) 


Now, 


Iff^GaVG-a)  -  E  llxllfp(u>*)  ~  E  II?* I! »*«) > 

*67(Gi)\7(G_i)  *67(Gi)\7(G_i) 


and  therefore  from  (3.26)  we  get 


\\PX-v\\hi(g3)  <  C'/i||p||^(Gi)||xl|H1(G3\G_2)> 
which  is  the  desired  result.  □ 

We  remark  that  the  Axiom  A3  as  stated  in  this  paper  is  slightly  different 
than  the  Axiom  A3  given  in  [30] .  We  further  remark  that  Axiom  A3  in  [30]  has 
been  used  only  to  prove  Lemma  2.3  (Page  911)  in  that  paper.  We  now  prove 
Lemma  2.3  of  [30]  using  the  Axiom  A3  as  stated  in  this  paper.  The  proof  is 
similar  to  the  proof  of  Proposition  2.2  in  [29]. 


Lemma  3.2  Let  D\  CC  D2  CC  D3.  There  exists  a  constant  C  such  that  given 
X  G  Sh(D3),  there  exists  rj  G  Sh(D3)  with  77  =  x  on  P2  such  that 

llx  -  r)Wm(D3\D2)  <  G||x||ffi(£,3\_Dl) 

and 

IMU,(r>3)  <  C'IIxI|l,(£>3)>  for  1  <  q  <  00. 

Proof:  Let  D\  CC  D2  CC  D-2 1  CC  D22  CC  D23  CC  D3.  Consider  p  G 
C°°(D22)  such  that  p  =  1  on  D23.  Then  from  (3.13)  and  (3.12),  with  D3  =  G3, 
D23  =  G2,  D22  =  G 1,  L>21  =  Go,  D2  =  G_1,  and  Di  =  G_2,  there  exists 
77  G  Sh{D3)  with  77  =  x  on  D2  such  that 

\\PX~v\\h^d3)  <  C'^IIxI|h1(d3\Di))  (3.27) 

and 

IMI Lq(D3)  <  G||x||l,(d3)- 
Since  p  =  1  on  D2 1,  we  have 

11(1  -  p)x\\m(D3)  <  C'||xI|j?1(b3\d2i)  <  C'||xl|i?1(u3\Di)) 
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and  thus  using  the  triangle  inequality  and  (3.27),  we  have 

IIX  -  VWH1  (D3)  <  11(1  -  p)x\\h1(d3)  +  II PX  -  vWh^Ds)  <  C'||x||h1(D3\£)1) 
Finally,  since  rj  =  x  on  D2,  we  have 

llx  -  v\\m(D3\D2)  =  ||x  -  v\\hi(d3)  <  C,||xI|h1(d3\£'i) 

which  is  the  desired  result.  □ 

Axiom  A4.  Scaling.  Let  the  sets  G  in  Axiom  Al,  G2  in  Axiom  A2,  G3 
in  Axiom  A3  be  the  sphere  Bp  CC  12  of  radius  D  >  C^h  with  center  xo-  The 
linear  transformation  y  =  (x  —  x0 )/D  takes  BD  into  a  sphere  B  and  Sh(BD) 
into  a  new  function  space  S(B).  Then  S(B)  satisfies  Axioms  Al,  A2  and  A3 
with  h  replaced  by  h/D.  Furthermore,  the  constants  occurring  in  Axioms  Al, 
A2,  and  A3  remain  unchanged,  in  particular,  independent  of  D. 

Using  a  standard  scaling  argument  argument,  one  can  show  that  Axiom  A4 
holds  with  respect  to  Axioms  Al  and  A2.  To  show  that  Axiom  A4  holds  with 
respect  to  Axiom  A3,  one  has  to  go  through  the  proof  of  Axiom  A3  with  Sh(BB) 
replaced  by  S(B).  We  skip  this  proof  in  this  paper. 

Axiom  A5.  There  exists  a  constant  C5  such  that  the  following  holds: 

(i)  For  any  Xq  G  such  that  the  ball  B 0  of  radius  h  centered  at  xo  is 
contained  in  f2,  there  exists  a  function  <5o  G  C 1  with  support  in  Bq  satisfying 

Xfao)  =  [  X?o,  for  all  x  e  Sh, 

Jb0 

and 

INU,  <  G5/i_Ar(1_1/9),  ||V50||l,  <  CBh-N^-1/q)-\  for  1  <  q  <  00.  (3.28) 

(ii)  Similarly,  for  j  =  1,  2, . . . ,  TV,  there  exists  5±j  such  that 

pL{x0)=  f  §^S-Lj  ,  for  all 

dXj  JBo  dx.j 

and  (3.28)  holds  with  So  replaced  by  Sij. 

We  will  show  the  existence  Si  j.  The  existence  of  <5o  can  be  shown  similarly. 
We  let  h  =  1.  Let  4)  be  a  smooth  non-negative  weight  function  with  compact 
support  in  Bq  and  suppose 


Consider  the  inner  product 


<  v,  w  >= 


dx. 
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and  let  il>i,ip2, . . .  ,il>e  be  an  orthonormal  basis,  with  respect  to  this  inner  prod¬ 
uct,  for  the  finite  dimensional  space  {J^-  :  %  G  Sh(B0)}.  Define 


e 

6l,j  = 

3=1 

Let  x£  Sh  such  that  ^  =  Em=i  cmVw  Then 

P  P 

/  ,jdx  =  y^lpjix o)  /  [V  Cm  1pm] Ipj $  dx 

Jb0  ax3  J=1  JB0  m=1 

=  '52'l>j(x0)cj  =  ^(*0)- 

3= 1  J 

We  get  (3.28)  by  scaling. 

Remark  3.3  We  mention  that  another  Axiom  AO,  related  to  a  trace  inequality, 
was  used  in  [30].  It  was  relevant  there  since  the  functions  in  Sh,  considered 
in  [30],  were  piecewise  C2  and  globally  W^.  In  our  case,  the  functions  in 
Sh  =  SGFEM  are  globally  C 2  and  thus,  a  trace  inequality  like  Axiom  AO  in  [30] 
is  not  needed. 


4  Superconvergence 

In  this  section  we  will  present  the  main  result  of  this  paper,  i.e.,  the  natural 
superconvergence  of  the  derivatives  of  the  Generalized  Finite  Element  solution 
in  the  interior  of  the  domain  12 ,  away  from  the  boundary  of  12.  The  analysis 
presented  in  this  section  will  closely  follow  the  analysis  given  in  [10],  [35]  in  the 
context  of  finite  element  method.  We  will  need  several  other  assumptions  that 
will  be  stated  in  this  section. 

Without  loss  of  generality,  we  assume  that  xq  =  (0,  0)  G  f2  and 

SA0  =  {a:  =  (x1,x2)  G  12  :  H^Hoo  =  max(|a;i|,  \x2\)  <2 H)  CC  12,  (4.1) 

where  H  >  0  will  be  determined  later.  We  also  define  the  set 

f2i  =  {x  G  12  :  Moo  <  H}.  (4.2) 

It  is  clear  that  the  solution  Uh  =  uqfem  of  (2.12)  satisfies 

B(u  -  uh,x)  =  0,  for  all  x  G  Sh(Q0),  (4.3) 

where  we  recall  that  Sh(Q o)  denotes  the  restrictions  of  the  functions  in  Sh(i 2) 
with  compact  support  in  the  interior  of  12q. 
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We  now  describe  several  assumptions  in  SCI  SC4  that  will  be  used  in 
the  analysis  presented  in  this  section. 

SCI:  We  assume  that  the  patches  u>i,  such  that  u;*  fl  fio  yf  0,  are  uniform 
and  translation  invariant  as  defined  below.  Consider  h  >  0,  where  h  is  of  the 
same  order  as  h  (recall  that  diam(wj)  <  2 h),  i.e. ,  there  exists  a  positive  constant 
C  such  that 

h  =  Ch. 

Let  i  =  (*1,^2),  where  ti,*2  =  0, ±1,±2, •••,  be  the  integer  lattice.  Suppose 
Xi  =  ih=  (iih,i2h)  is  the  center  of  the  circle  lo*  C  u>i  for  i  £  7(^0),  where 

7(f2o)  =  {i  £  Z  x  Z  :  w.j  n  Oo  7^  0} 

We  note  that  we  are  enumerating  u>i  differently  from  the  way  we  enumerated 
them  in  (2.5).  Clearly,  Xi+j  =  Xi  +  Xj  for  i,j,  i+  j  £  7(^0)-  We  further  assume 


that  for  all 

i  £  7(^0), 

uii  =  {x  +  Xi  :  x  £  wo}. 

(4.4) 

uj*  =  {x  +  Xi  :  x  £  tv  J}, 

(4.5) 

4>i(x)  =  4>o(x  -  Xi),  X  £  U)it 

(4.6) 

Vi  =  |fi(®)  :  &(%)  —  £o(x  —  Xi),  where  £0  €  Vo}. 

(4.7) 

Clearly, 

<j>i(x  —  Xj)  =  (t>i+j{x),  for  i,j,  i  +  j  £  7(^0) 

(4.8) 

and 

Vi+j  =  {&(*  -  Xj)  :  £i  £  Vi },  for  i,j,i  +  j  £  7(^0) 

(4.9) 

SC2:  We  consider  H  in  (4.1)  and  (4.2)  of  the  form 

II 

(4.10) 

where  0  <  f3  <  1  will  be  determined  later.  We  let 

M0  =  {i£ll:  a;  |oo  <  h/2) 

(4.11) 

and  define 

Mj  =  {x  £  ft  :  x  =  Xj  +  y,  where  y  €  M0} 

(4.12) 

Clearly, 

Mj  =  {x  £  0  :  || ar  —  £j||oo  <  /i/2} 

(4.13) 

and  we  assume  the  h  and  H  are  such  that 

Ll0=  U  Mj. 

(4.14) 

XjGfio 

SC3:  For  a  given  H ,  let  H  >  0  such  that 

Ho  —  |  J  (jJj  CZ 

je-y^o) 
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where 


S]o  =  {^fi:  Moo  <H}  (4.15) 

and  H  <  CH  where  1^<  C  is  a  fixed  constant.  We  assume  that  the  solution  u 
of  (2.4)  is  smooth  in  Qq,  i.e., 

S  c  <446> 

where  k  is  the  degree  of  the  polynomials  in  V)  =  Vk(uij). 

SC 4:  We  assume  that  the  solutions  u  and  Uh  of  (2.4)  and  (2.12)  satisfy 

Chk  <  \u  —  Uhlw^n-t)  (4-17) 

\\u-uh\\LooiQo)<Chk+1-\  (4.18) 

where  0  <  e  <  1. 

We  now  define  a  function  ph(x),  which  will  play  an  important  role  in  the 
analysis  presented  in  this  section.  Towards  this  end,  we  first  define  a  linear 
operator  ij*  :  W^+1(u;j)  — >  V)  =  Vk(u Ji)  for  each  i  G  7(^0)  satisfying 

( i )  Ii\pk(-)](x)  =  Pk(x),  for  x  G  tot,  and  for  all  pk  G  Vk  (4.19) 

(m)  For  v  G  W'(^+1(w,:), 

II*  -  <  cdk+1-l\\v\\w^{UiV 

for  all  0  <  l  <  k  +  1  and  1  <  q  <  00  (4.20) 

(in)  For  v  G  W£+1(bJj+i)  and  for  all  1,1+  j  G  7(flo) 

Ik+j[v(-)](x  +  Xj)  =  Ii  [v(-  +  Xj))(x),  for  all  x  G  W|.  (4.21) 

For  /f[u(-)](a;),  one  could  take  the  restriction  of  the  Taylor  polynomial  of  v(x), 
centered  at  Xi,  to  u>i-  We  then  define  the  operator  Ih  :  W^+1(f2o)  — 1 ►  Sh(Slo)  by 

ih[Q(-)}(x)=  M*)ii[Q{-)}(x),  xen0,  (4.22) 

*6  7(f20) 

where  Q  G  W'(^+1(f2o)-  Clearly,  for  a  polynomial  pk  of  degree  k, 


ih\pu(-)](x) 

=  £  Mx)Ii\Pk(-)\(x) 

iS  7(f2o) 

=  ^2  4>i(x)pk(x)  =  pk(x),  x€fl0  (4.23) 

ie7(fio) 

The  operator  Ih  also  satisfies  the  standard  interpolation  estimate  given  in 
the  following  lemma. 
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Lemma  4.1  Let  v  G  W^+1(flo)-  Then 

ll'W-^K)]  Ilw*(n0)  <  Chk+1~tWvWw^(floy  forl<q<oo  and  t  =  0,1,2. 

(4.24) 


Proof:  We  first  note  that  for  x  £  LIq, 

v-ihv  =  v-  Y  & 4X-)]  =  Y  Mv - Jlib(-)])- 

*67(^0)  *67(^0) 


Therefore,  using  (2.6),  (2.9),  and  (4.20),  we  have 


v.)  s  c  e  t*(»--r?w-)i)r„.(„,) 

ie^C^o) 

t 


<  c  y.  Ew 


*6 7(°o)  3—° 


Q 

wL(uh) 


v-If\v 


Wl~3{ui 


<  c  y  < 

*£7(^0) 

<  Ch(k+1-t)q\\v\\q 


(k+l  —  t)qn  nq 


'Wf  +  1(uH 


Wf+\Q.  0) 


where  we  used  the  fact  that  di  <  2 h  in  the  last  step.  Thus  we  get,  for  small  h, 

h-iX-Mwu^^ch^M 


'wf+1(n0y 


which  is  the  desired  result.  □ 

We  finally  define  ph(x)  as 

ph(x)  =  Q(x)  -  Ih[Q(-)](x),  x  G  fl0,  (4.25) 

where  Q{x)  is  a  polynomial  of  degree  k  +  1. 


Lemma  4.2  ph(x)  is  periodic  in  f2o>  *.e., 

ph(x)  =  ph{x  +  Xj),  for  x  G  M0  and  x  +  Xj  G  O0. 

Proof:  We  first  note  that 

Q(x)  =  Q(x  +  Xj)  —  pk(x;xj),  (4.26) 

where  pk{x\  Xj)  is  a  polynomial  of  degree  k  that  depends  on  Xj.  Now  from  the 
definition  of  ph(x)  and  using  (4.26),  (4.23),  we  have 

A  x)  =  Q{x)  -  Ih[Q{-)]{x) 

=  Q(x  +  Xj)  -  pk(x;  Xj)  -  Ih[Q(-)](x) 

=  Q(x  +  Xj)  -  Ih\pk(-,Xj)  +  <3(-)](®) 

=  Q(x  +  Xj)  -  Ih  [Q{-  +  Xj)](x).  (4.27) 
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We  will  now  show  that 


Ih[Q( ■  +  Xj)](x)  =  Ih[Q(-)]( x  +  Xj),  for  x  £  Mo  and  x  +  Xj  €  fio- 
Since  x  £  M0,  we  have  x  +  Xj  £  Mj,  and  using  (4.21),  we  get 

Ih[Q(-)](x  +  Xj)  =  Y  (t>i{x  +  x3)Ii[Q{-)}{x  +  xj) 

iei(Mj) 

Y  +  xi) 

iei(Mj) 

Y  <t>i(x)tf+AQ(')  }(x  +  xi) 

ie~r(M0) 

=  Y  +  xj)](x) 

lej(Mo) 

Thus  from  (4.27)  we  get 

ph{x)  =  Q(x  +  Xj)  -  Ih[Q{-)]{x  +  Xj)  =  ph(x  +  Xj) 
which  is  the  desired  result.  □ 

The  above  result  can  also  be  stated  as  ph(x)  £  Hpol.(n0),  where 

Hper(fi0)  =  {»£  B1(fl0)  nC°(fi0)  :  v(x)  =  v(x  +  Xj )  for  x  £  AI0lx  +  Xj  £  f20} - 

(4.28) 

Next,  for  a  given  v  £  Hpel.(idp),  we  will  define  its  periodic  projection 
PpeiV  £  Sh(n0)  fl  Hper(tt0),  and  present  two  results  that  will  also  be  used 
in  our  analysis.  We  first  consider  the  subspace 

Spfter(M0)  =  {xe  sh  (M0)  :  x(-h/2,  x2)  =  X(h/ 2,  x2), 

x(xi,-h/2)  =  x{xi,h/2), 

for  |aii|  <  h/2,  |cc2 1  <  h/ 2}.  (4.29) 

For  v  £  Hper(fl0),  we  define  Pperv  £  Sper(M0)  as  the  projection 

BMo(PperV,X)  =  Bm0  ( v,x ),  for  all  x  £  S%ei(M0),  (4.30) 

[  (v  -  PpeIv)  dx  =  0,  (4-31) 

J  Mo 

where  Bm0{u,v)  =  fMo  Vu-Vvdx  (see  2.3)).  We  then  extend  Pperv  periodically 
to  ft0,  again  denoted  by  Pperv,  as 

Ppelv(x)  =  Pperv(x  —  Xj),  for  all  x  £  fio,  x  —  Xj  £  Mq-  (4.32) 

Thus  Ppeiv  £  Sh(flo)  fl  Hpei(Q0). 
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(4.33) 


Lemma  4.3  Let  v  £  HpCr(flo).  Then 

B(v  -  PPerV,  x)  =  0,  for  all  \  G  Sh(Cl 0). 

Proof:  Since  Mj’s  (i.e.,  interior  of  Mj)  are  non-intersecting  and 
flo  =  Ux.en 0Mj,  we  note  that 

B{v  -  Pper^,  x)  =  5Z  BmJ  “  -Pper^,  X)  •  (4-34) 

Now  using  the  periodicity  of  v  —  Pperv,  we  get 

BMj  (V  -  PperV,  X)  =  /  V(U  -  Pper^)  '  Vx  dx 

Jm, 

=  /  v(w(j/  +  Sj)  -  Ppeiv(y  +  Xj))  ■  Vx{y  +  Xj)  dy 

J  M0 

=  [  V (u(y)  -  Pper^(y))  •  Vx(y  +  arj)  dy  (4.35) 

J  M0 

Thus  from  (4.34),  we  get 

-  ^per^,x)  =  /  v(v  -  Pper^)  -^xdy  =  BMq(v  -  Pperv,x)  (4.36) 
J  Mo 


where 


x(y)=  x(y  +  Xj). 


Since  x  G  Sh(Llo),  we  can  show  that  %  £  Sper(Mo).  Thus  using  (4.30)  in  (4.36), 
we  get 


B(v  -  PperV,  x)  =  o 


which  is  the  desired  result.  □ 


Lemma  4.4  Let  v  £  Hper(flo).  Then 

||V  -  PperV|U2(fi0)  <  Ch\\v\\Hi(Q0y  (4.37 

Proof:  We  first  note  that  since  v  —  Pperu  is  periodic,  using  (4.31)  we  have 

/  (v-  Pperv)  dx=  (v-  Pperv)  dy  =  0, 

J  Mj  J  Mo 

and  using  Poincare  inequality,  we  get 

||^  Pper^  ||  L2  {Mj  )  ^  Ch\V  Pper'^|ff1(MJ)  h\v  Pper'^|//1(Mj)  • 
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Therefore, 


11^  —  PperV||i2(n0)  —  ^  llW  “  -PperV||  L2(Mj) 

Xj£Q  o 

—  Ch,  ^  K;  —  -Pper^l  H1(Mj) 

XjEflo 

-  Ch  ^2  +  \PpeM2Hi(Mj)} 

XjEfi,0 

=  Ch{\v\2Hi(ao}  +  |-PPer^|l/i(a0)}-  (4.38) 

Again  using  the  periodicity  of  v  and  Ppelv,  and  the  fact  that  Ppelv  |M  is  the 
projection  of  v\Mq  onto  S^er(AI0),  we  have 

| PperV |  H1  (Oo)  <  C|  v\h'  (fi0) » 

and  thus  from  (4.38),  we  get  the  desired  result.  □ 

We  now  present  our  main  theorem. 

Theorem  4.1  Suppose  the  assumptions  A1  A3  hold  and  the  assumption  SCI 
SC4  are  satisfied  with  0  <  /3  <  1  —  e,  where  j3  and  e  are  as  in  (4-10)  and  (4-18) 
respectively.  Also  assume  that  X^|s|=fc+i  l^s'u(a:o)|2  >  0,  where  s  =  (si,S2)  is  a 
multi-index  and  Xq  =  (0,0).  Then,  for  h  small  enough,  there  exists  a  >  0  such 
that,  for  i  =  1,2, 

—  (u- uh)(x)  =  —  (ph  -  Pperph)(x)  +  Ri(x),  forxefli,  (4.39) 
where 

lli^wno  <  C7»fc+“, 

and  Q(x)  in  the  definition  of  ph(x)  (see  (4-25))  is  the  (k  +  l)th  degree  Taylor 
polynomial  of  u  centered  at  Xq  ■ 

Remark  4.1  We  note  that  we  have  assumed  u  to  be  smooth;  in  particular,  u 
satisfies  (4.16)  in  assumption  SC4- 

Remark  4.2  If  a;*  €  fii  is  a  zero  of  (ph  —  Pperph),  then  from  the  above 
result  we  get 

J^-(u-uh)(x*)  <  Chk+a, 
and  thus  x*  is  a  natural  superconvergence  point  of  ( u  —  Uh). 

Proof  of  Theorem  4-1:  The  proof  will  be  given  in  four  steps. 
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1.  Since  Q  is  the  (k  +  l)th  degree  Taylor  polynomial  of  u  centered  at  Xq  = 
(0, 0),  we  have 

llu“  Qllw^tfio)  <  \\u  -  QWw,  (^o) 

<  CHk+2~s  <  CHk+2~s,  for  0  <  s  <  fc  +  2,  (4.40) 

where  Q0  is  defined  in  (4.15).  We  now  define  the  Neumann  projection  PfcQ  of 
Q  onto  Sh((l0)  by 

Bno(Q-P%Q,x)  =  0,  for  all  X  €  Sh(Si0)  (4.41) 

[  (Q-PkQ)dx  =  0.  (4.42) 

We  write 

U~Uh  =  ( Q  —  PjyQ)  +  [(ti  —  Q)  —  ( llh  —  PjyQ)^  , 

and  thus, 

■S—  (u  -  Uh)(x)  =  (Q  —  PmQ)(x)  +  n(x),  for  x  £  fii,  (4.43) 

OXi  OXi 

where 

IkilUooCni)  <  I (u  —  Q)  —  (uh~  PnQ)\w^{Q!)-  (4.44) 

We  will  estimate  the  term  on  right  hand  side  of  this  inequality. 

2.  Since  Sh(flo)  C  Sh(Qo),  from  (4.3)  and  (4.41)  we  have 

B((u- Q)  -  (uh  -  PffQ),x)  =  0,  for  all  x  e  Sh(ClQ).  (4.45) 

Therefore  using  (4.44)  together  with  the  interior  estimate  (3.2)  with  D  =  H, 
Qd  =  Clo,  and  Qy  in  (3.2)  replaced  by  fii  in  (4.2),  we  have 

IkilUoc^!)  <  C  min  [||(u-Q)-x||lvij(no)  +  H'1||(u-Q)-x||ioo(no) 
xesh(Q o)  L  J 

+CH-2\\(u  -Q)-  (uh  -  PnQ)\\l2(q0)-  (4-46) 

From  the  approximation  axiom  (3.3)  and  (4.40),  we  get 

ll(«  -Q)~  Xlki(no)  <  Chk\\u  -  Q\\w^o)  <  ChkH. 

Similarly,  we  get 

\\(u-Q)-x\\L^0)<Chk+1H. 

Therefore,  from  (4.46)  we  have 

INIloo(Qi)  <  ChkH  +  CH-2\\{u  -  Q)  -  (uh  -  P*Q) |U2(no),  (4.47) 

where  we  used  the  fact  that  h  <  H. 
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Now  using  the  assumption  (4.18),  we  get 


Il«  -  Uhh2(a0)  <  CH\\u  -  WfcllL^n,,)  <  CHhk+1  £.  (4.48) 

Since  Hq  is  convex,  using  a  standard  duality  argument,  a  property  of  the  pro¬ 
jection  P^Q,  and  (4.40),  we  get 

ii<5-cteiu,(n„)  <  ch\Q  - 

<  ch^Hmw^&) 

<  C/.i+,H||»||Hi„(Bii)  <  Chw-H.  (4.49) 

Therefore,  from  (4.47)  and  (4.48)  we  have 

lki||i00(n1)  <  ChkH  +  CH  2  ||w  —  Uh||i2(n0)  +  ||Q  —  P^Q||L2(n0) 

<  ChkH  +  CH~1hk+1~e.  (4.50) 

3.  We  now  consider  the  term  ( Q  —  P^Q)  in  (4.43).  Let 

i/j(x)  =  ph{x )  -  Pperph{x), 

where  ph( x)  was  defined  in  (4.25)  (with  Q  as  in  this  theorem)  and  PpeTph(x)  is 
its  periodic  projection  (see  (4.30)— (4.31)) .  We  write 

Q  -  PnQ  =  4>  +  Q-  PnQ  -  i/>. 

Thus,  for  x  €  fir, 

d  d 

—  (Q  -  PxQ)(x)  = —ip(x) +  ri(x),  (4.51) 

where 

||f:i||L00(n1)  <  | Q  —  PnQ  ~  V’lwicni)-  (4-52) 

To  estimate  the  right  hand  side  of  the  above  inequality,  we  note  that 

Q  —  PnQ  ~  4’  =  Q  ~  PnQ  ~  ph  +  Ppei-ph 

=  Q-PkQ-Q+Ih{Q(-)}+PperPh 
=  Ih[Q(-)]-PkQ  +  PpeiPh, 

and  thus  Q  —  P^Q  —  ip  €  Sh(fl0). 

We  recall  that  ph  €  L7per(flo).  Therefore  from  (4.33)  and  (4.41),  we  get 
B(Q-PkQ-iP,X)=  0,  for  all  %  €  Sh(flo). 
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Hence  from  the  interior  estimate  (3.2)  with  u  =  0  (also  D ,  fir),  and  Ho  redefined 
as  before)  and  using  (4.49)  and  (4.37),  we  get 


\Q  -  PnQ  -  V’lv^cni)  ^  c#  2\\Q  -  pnQ  -  V’IImho) 

<  Cff-2||Q  -  PnQ\\l2(q0)  +  C'H-2||^|U2(n0) 

<  CH-'h^1-*  +CH-2h\\ph\\HHno).  (4.53) 


Also  from  (4.24)  and  (4.40),  we  have 

ll/^lliTH^o)  =  IIQ  -  Ih[Q(')]  II^HOo) 
<  Chk\\Ql 


fffc+1(0 o) 


<  ChkH\\Q\\ 


wd+1(n0) 


<  ChkH\\u\\ 


'WSo+1(Q0) 


<  ChkH. 


Thus  from  (4.52)  and  (4.53),  we  get 

INUoo(ffi)  <  I Q-  PnQ  ~  V’lwitno  <  CH~1hk+1~e  +  CH~lhk+1.  (4.54) 

4.  Finally  combining  (4.43),  (4.50),  (4.51),  and  (4.54),  and  writing  Rpx )  = 
r*(x)  +  fi{x),  we  obtain 

d  d 

——  (u  -  uh){x)  =  ip(x)  +  Ri(x) j  for  x  G  Hi, 
ax*  oxi 

where 


ll-RilUooCOO  <  IkilUcoino  +  INUoctfiO  <  ChkH  +  CH~1hk+1~e  +  CH-1hk+1. 

We  now  recall  that  H  =  hP ,  where  h  =  Ch  and  0  <  (3  <  1  to  be  determined. 
Therefore, 

llifclli^no  <  Chk+0  +  Chk+1-0-e  <  +  h 1~0~e) 

We  choose  (3  such  that  0  <  /3  <  1  —  e,  and  define 


a  =  min  (/?,  1  —  j3  —  e), 


to  get 

llifcllwno  <  Chk+a, 

which  is  the  desired  result.  □ 


Remark  4.3  In  the  proof  of  Theorem  4.1,  we  considered  ph(x)  =  Q{x)  — 
Ih[Q(-)](x),  where  Q(x)  was  the  Taylor  polynomial  of  degree  k  +  1  centered 
at  xq.  In  fact,  it  can  be  easily  shown  that  (4.39)  holds  when  Q(x)  is  the  poly¬ 
nomial  ^f=o1[^i^2+1_*u(0i  0)]  XiX2+1_*  of  degree  k  +  1  (a  linear  combination  of 
monomials  of  degree  (k  +  1)). 
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Remark  4.4  In  Remark  4.2,  we  have  seen  that  the  zeros  of  (ph—Ppelph)  are 
the  superconvergence  points.  Because  of  the  periodicity  of  ph  —  PpeTph,  we  only 
need  the  zeros  xjj  of  ^  (ph  —  Pperph)  in  the  cell  M0.  All  other  superconvergence 
points  in  Q 1  can  be  found  by  a  simple  translation  x*  =  xp+Xi,  where  xp  +  Xi  G 
Q-i .  The  points  Xq  can  be  obtained  by  finding  the  zeros  Xq  of  g|-  ( ph  —  Pperph) 
in  the  “master  cell”  with  h  =  1,  and  then  scaling  back  xp  to  Mq.  Thus  the 
computation  of  x*  does  not  depend  on  h  or  the  solution  u  of  (2.4). 


Remark  4.5  It  is  immediate  from  (4.17)  and  (4.39)  that 

d  r 


and 


dx  Uh ^  ^  PperP  )J  ||^oo(f21)  —  ^  \U  Uh\wl0{Sli) 

d 


dxi 


( u  -  uh)(x*)  <  ha[u  - 


where  x*  is  a  superconvergence  point. 


We  recall  that  we  considered  the  partition  of  unity  functions  <f>i,  used  in 
gGFEM ,  p,e  (j2  fLinctions,  and  thus  Uh  G  c/GFEM  jg  aiso  a  pj2  function. 
We  will  now  present  a  result  on  the  superconvergence  related  to  the  second 
derivatives  of  u  —  Uh ■  The  analysis  will  be  essentially  same  as  the  analysis  in 
Theorem  4.1,  but  we  will  require  additional  assumptions. 

Let  H2  C  Hi  be  a  square  centered  at  Xo  given  by 

H2  =  {xGH:  ||a:||oo  <H/2}. 

In  addition  to  the  assumptions  in  Theorem  4.1,  we  assume  that  k  >  2  and 

Chk~1  <  \\u  —  Uh\\w^(n2)  (4.55) 

We  now  present  a  theorem,  which  is  another  important  result  of  this  section. 
In  the  proof  of  this  theorem,  we  will  use  certain  technical  results  obtained  in 
the  proof  of  Theorem  4.1. 

Theorem  4.2  Suppose  all  the  assumptions  of  Theorem  f.l  hold.  Also  suppose 
that  (4-55)  is  satisfied  and  k  >  2.  Let  s  =  (si,s2)  be  the  multi-index  with 
|s|  =  2.  Then ,  for  h  small  enough,  there  exists  a  >  0  such  that 

Ds(u  -  uh)(x)  =  Ds(ph  -  Ppeiph)(x)  +  Rs(x),  x  G  fl2,  (4.56) 


where 

IlflJwn,)  <  Chk~1+a  <  Cha\\u-uh\\wi(n2). 

and  Q(x)  in  the  definition  of  ph(x)  (see  (4-25))  is  the  (k  +  l)th  degree  Taylor 
polynomial  of  u  centered  at  Xq  =  (0,0). 
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Remark  4.6  We  note  that,  as  in  Theorem  4.1,  we  have  also  assumed  u  to  be 
smooth  in  this  theorem  ;  in  particular,  u  satisfies  (4.16)  in  assumption  SC4- 


Proof:  With  Q(x)  and  P'fQ  as  in  the  proof  of  Theorem  4.1,  we  write 

Ds(u  -  Uh)(x)  =  DS(Q  -  PjfQ)(x)  +  rs(x),  for  x  e  ^2,  (4-57) 

where 

llrslUoo(o2)  <  || («  -  Q)  -  ( Uh  -  PnQ)\\w^(Q2)-  (4.58) 

Let  Vh  G  Sh(flo)  be  arbitrary.  Using  the  inverse  inequality  (3.7),  we  have 

|| (u  —  Q)  -  (uh  —  PnQ) ||w£,(fi2)  <  \\iu  ^  Q)  ~  vh\\w^(n2) 

+  lluft  _  iuh  ~  -f>w<3)llwr^(n2) 

<  \\{u  -  Q)  ~  Vh\\w^(n2) 

+Ch  1\\vh  —  (uh  —  P^Q)  llw^(Qi) 

<  \\{u  -  Q)  ~  Vh\\w^(n2) 

+Ch  1 1|  (it  —  Q)  —  v/tllw^ni) 

+Ch  1\\(u  -  Q)  —  (uh  -  PnQ) ||w^,(ni)- 

Therefore,  using  the  approximation  property  (3.3)  and  (4.40)  in  the  proof  of 
Theorem  4.1,  we  get 

\\(u-Q)~(uh-  PhQ)\\wl(n2)  <  Chk~l\\u- QWyyk+i^ 

+Ch~1\\(u  -  Q)  -  ( Uh  -  -P^Q)Hrv^(ni) 

<  CHh fc_1 

+Ch~1\\(u  -  Q)  -  ( Uh  -  -P^Q)||rv(jo(n1)- 

(4.59) 


A  careful  examination  of  the  arguments  leading  to  (4.50)  shows  that 
||(u  -Q)~  ( uh  -  ^<2)11^1,(00  <  ChkH  + 
and  thus  from  (4.58)  and  (4.59),  we  get 

IN|Loo(n2)  <  CHh*-1  +  CH-1hk~e.  (4.60) 

Again,  with  xf>  as  in  the  proof  of  Theorem  4.1,  we  write 

DS(Q  -  PnQ)(x)  =  Dsif  +  rs,  for  x  G  Cl2,  (4-61) 

where 

ll^llLoo^a)  <  IIQ  -  pnQ  -  V’ll w£,(o2)-  (4.62) 

We  have  seen  in  the  proof  of  Theorem  4.1  (in  the  paragraph  after  (4.52))  that 
Q  —  P^Q  —  4>  &  Sh(fl 0).  Therefore  using  the  inverse  inequality  (3.7)  we  get, 

||Q  —  PnQ  ~  ^llw^(a2)  <  Ch  1 1| Q  -  PnQ  —  "011  (sii) •  (4.63) 
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We  have  also  shown  in  (4.54)  in  the  proof  of  Theorem  4.1  that 


II Q  -  PnQ  -  n  w^n,)  <  CH~1hk+1^e  +  CH~1hk+1, 

and  thus  from  (4.62)  and  (4.63),  we  get 

||f8||ioo(n2)  <  CH~xHk~*  +  CH~lhk.  (4.64) 

Finally,  writing  Rs(x)  =  rs(x)  +  ra(x)  and  combining  (4.57),  (4.60),  (4.61), 
and  (4.64),  we  get 

Ds(u  —  Uh)(x )  =  Dstp(x )  +  Rs(x),  for  x  &  0.2, 

where 

||-RS|Uoo(f22)  ^  IkslUoctOa)  +  ||f;s|U00(Q2) 

<  CHhk~l  +  CH~1hk~c 
=  Ch^^h0  +  h1~0~e) 

Thus,  choosing  0  <  (3  <  1  —  e  and  a  =  min (/?,  1  —  /?  —  e),  and  using  (4.55),  we 
get 

IlflJwno  <  Chk~1+a  <  Cha\\u-uh\\wl{n2), 
which  is  the  desired  result.  □ 


Remark  4.7  We  note  that  following  the  arguments  presented  in  the  proof  of 
Theorem  4.2,  it  is  possible  to  obtain  a  superconvergence  result  like  (4.56)  for 
higher  derivatives  of  u  —  Uh,  i.e.,  for  Ds(u  —  Uh)  for  |s|  >  2.  We  do  not  give  a 
proof  this  result  here  to  keep  the  exposition  simpler. 


5  Example 


In  this  section,  we  will  present  a  computational  example  to  illuminate  the  results 
given  in  Section  4. 

We  consider  the  one-dimensional  version  of  the  problem  (2.1)  with  O  =  (0, 1), 
where  the  exact  solution  is  u(x)  =  sin(7ra;/2).  To  approximate  this  solution  by 
GFEM,  we  choose  nodes  Xi  =  ih,  i  =  0, 1,  •  •  •  ,  N,  where  Nh  =  1,  and  define 
patches  uy  as  in  Example  1  in  Section  2  (see  (2.25),  (2.24)).  For  partition  of 
unity  functions  to  be  used  in  the  GFEM,  we  employ  functions  <f>i(x),  as  defined 
in  (2.27),  with  r  =  0.3  and  s(x)  =  q(x  —  r),  where 


q(y) 


0  <  y  <  1  —  2r. 


We  use  the  space  of  linear  polynomials  for  local  approximating  spaces,  i.e.,  Vj  = 
7?1(wj).  We  denote  the  GFEM  approximation  of  u(x)  by  ugfem{x)  =  Uh(x). 
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It  is  clear  from  (4.39)  of  Theorem  4.1  that,  for  h  small,  the  superconvergence 
points  of  v!  —  u'h  in  [xj,Xj+ 1]  CC  O  are  the  roots  of  \ph  —  PpeTph]'  in  [xj,xj+ 1]. 
We  obtain  these  roots  by  first  finding  the  roots  y*  of  ^  [p—  Ppelp]  on  the  master 
“cell”  0  <  y  <  1.  Here 

l 

p  =  y2 

i= 0 

where  <j>i(y)  is  the  PU  function  with  h  =  1  and  Ii(y)  is  the  Taylor  polynomial  of 
y2,  restricted  to  w,;  with  h  =  1.  Also  Ppelp  is  defined  as  Pperp  €  <SPer,  such  that 


l 


Pper/S  dy 


for  all  V  £  S'per 


where 

S'per  =  spanji,  4>o{y)y  +  tyiidijiy  -  !)} 

Finally,  the  superconvergence  points  x*  of  v!  —  v!h  in  [xj,  Xj+h]  is  given  by  scaling 
as 

x*  =  Xj  +  y*h. 


Figure  5.1:  Graph  of  p  —  Pperp]  on  the  master  cell  [0, 1]. 

In  Figure  5.1,  we  present  the  graph  of  j^[p  —  Pperp]  in  [0, 1].  The  roots  of 
this  function  are  y{  =  0.058309  and  y%  =  0.555216.  Consequently,  the  super¬ 
convergence  points  of  v!  —  v!h  in  [xj,Xj+ 1]  are 

x*  =  Xj  +  0.058309  h,  and  x*2  =  x:j  +  0.555216  h.  (5.1) 
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(a)  (b) 

Figures  5.2  (a)  and  (b):  (a)  Graph  of  u'  —  u'h  on  Q  =  (0, 1)  with  h  =  0.1. 

(b)  Graph  of  u'  —  u'h  on  (xj,Xj+ 1)  =  (0.5,  0.6) 

In  Figure  5.2a,  we  present  the  graph  of  the  error  v!  —  u'h  on  O  =  (0,1), 
where  Uh  is  the  GFEM  approximation  of  u  with  h  =  0.1.  It  is  interesting  to 
note  that  u'  —  u'h  is  zero  at  several  points  in  the  domain  O  =  (0, 1).  In  Figure 
5.2b,  we  present  the  graph  of  u'  —  u'h  on  (xj,Xj+ 1)  =  (0.5, 0.6),  and  also  show 
the  superconvergence  points  x\  =  0.5058309  and  x%  =  0.5555216.  It  is  clear 
from  Figure  5.2b  that  |(m'  —  u'h)(x*)\  for  i  =  1,2  is  much  smaller  than  the 
max  [(«'  —  u'h)(x)  |,  .5  <  x  <  .6. 

We  next  computed  the  GFEM  approximation  Uh  for  h  =  0.1,0.05,0.025, 
and  0.0125.  For  each  value  of  h,  in  Table  5.1  we  display  M  =  ma x(u'  —  u'h)(x), 
x  £  [xj,Xj+ 1]  =  [0.5, 0.5  +  h],  e'  =  | (v!  —  ^4)(;ci)l  anc^  for  i  =  1,2,  where 
x*  are  the  superconvergence  points  in  [0.5, 0.5  +  h],  given  in  (5.1). 


h 

M 

ei 

e'i/M 

e2 

e'2/M 

0.1 

1.13  x  10"1 

2.68  x  lO"3 

2.37  x  10“2 

1.89  x  10"3 

1.66  x  10~2 

0.05 

5.42  x  10"2 

6.76  x  10-4 

1.25  x  10"2 

4.97  x  10"4 

9.16  x  10"3 

0.025 

2.65  x  10"2 

1.70  x  10"4 

6.41  x  10"3 

1.27  x  10"4 

4.80  x  10"3 

0.0125 

1.31  x  10"2 

4.30  x  10"5 

3.29  x  10"3 

3.15  x  10"5 

2.41  x  10"3 

Table  5.1 


It  is  clear  from  Table  5.1  that  the  ratios  e\ / M  and  e'2 / M  decrease  as  h  decreases, 
which  illuminates  the  Remark  4.5.  It  also  indicates  that  x\  and  a;2  are  indeed 
superconvergent  points  of  v!  —  u'h  in  [xj,Xj+ 1]  =  [0.5, 0.5  +  h]. 
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